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Abstract We analyze geometrical structures necessary to represent bulk 

and surface interactions of standard and substructural nature in complex 

bodies. Our attention is mainly focused on the influence of diffuse interfaces 

rl^ . on sharp discontinuity surfaces. In analyzing this phenomenon, we prove the 

^H| covariance of surface balances of standard and substructural interactions. 
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1 Introduction 

Bodies are called complex when their material substructure (i.e. the texture 
from nano-level to meso-level) has a prominent influence on their gross be- 
havior. Such an influence appears also through a not negligible occurrence of 
interactions due to substructural changes. Examples are liquid crystals, elas- 
tomers, ferroelectric and microcracked bodies, spin glasses. Above all, soft 
condensed matter displays complex behavior. Applications in nanotechnol- 
ogy, smart structures and various fields of technology open basic theoretical 
and experimental problems that challenge, in a certain sense, even some 
aspects of the foundational concepts of standard continuum mechanics. Ba- 
sically, the standard paradigm of Cauchy's format of continuum mechanics, 
prescribing that the material clement is a sort of indistinct sphere that we 
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collapse in a point in space seems to be not sufficient to account for the 
articulated substructural nature of a complex body. 

In fact, for complex bodies the material element is rather a 'system' and 
one needs the introduction of appropriate morphological descriptors u of 
such a system (order parameters) at least at a coarse grained level. They 
describe the essential geometrical features of substructural shapes. 

Physical circumstances of disparate nature suggest a rich crop of possible 
choices of v, each one characterizing special models. Moreover, the selection 
of morphological descriptors is strongly related with the representation of 
substructural interactions arising within each material element and between 
neighboring material elements as a consequence of substructural changes. 
Interactions are represented in fact by objects conjugated in the sense of 
power with the rates of the quantities describing the geometry of the body 
and its changes. In this sense, since placement and order parameter fields 
are involved, the description of complex bodies adopted here is called mul- 
tifield. It has basic differences with standard internal variable models. In 
a multifield approach, morphological descriptors enter directly the geomet- 
rical representation of the body and its kinematics; true interactions are 
associated with their rates and balanced. Information about the material 
substructure are introduced already at the level of geometrical description 
of the body. On the contrary, in standard internal variable models, the ge- 
ometrical description of the body is of Cauchy's type: the material element 
is morphologically equivalent to an indistinct sphere described just by its 
place in space. There, internal variables come into play in constructing ther- 
modynamics to carry information about the material substructure just at 
constitutive level (a detailed treatment of these classes of models can be 
found in [34]). 

By following the general unified framework of multifield theories pro- 
posed by Capriz in 1989 [4] (see also [3]) and then developed further on in 
its abstract structure (see [32], [23], [8], [9]), we do not specify the nature 
of v. We require only that v be an element of a finite-dimensional differen- 
tiable paracompact manifold without boundary M. to cover a class as large 
as possible of special theories. Our attention is then focused on the gen- 
eral setting which contains as special cases prominent theories interpreting 
problems typical of condensed matter. In a certain sense our work deals 
with a model of models. However, if we tackle this general point of view, 
we face the basic difficulty that M. does not coincide with a linear space 
in general. Moreover, we cannot consider a priori M. embedded in a linear 
space because the embedding is a non-trivial part of some special modeling. 

We focus our attention on conservative processes. For them, the relevant 
appropriate Hamiltonian formalism has been developed in [9] as an evolu- 
tion of the Hamiltonian formalism in classical non-linear elasticity (the one 
discussed in [27]). We start from the results in [9] and analyze some of the 
rather subtle geometrical questions induced by the abstract nature of M.. 

Our essential point of view is the one of [25] : " Geometry and mechanics 
of maps between manifolds are a general framework for condensed matter 
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physics and are also a tool to construct new models of unusual and perhaps 
unaspected phenomena." 

The main attention is focused on the representation of the way in which 
diffuse interfaces influence the behavior of additional sharp interfaces. In 
fact, the presence of the gradient of v in the list of entries of the energy 
may account for diffuse interfaces as well as it may describe weakly non- 
local (i.e. gradient) effects in homophase materials. Additional discontinuity 
surfaces may also occur as a consequence of the presence of defects, shock 
or acceleration waves, solidification phenomena of complex fluids etc. As 
an example one may consider a polarized ferroelectric material in which 
external loads induce a shock wave: the shock front encounters walls of 
polarized domains and interact with them. So, in general, one may say 
that, in principle, diffuse interfaces and potential homophase gradient effects 
may interact with sharp discontinuity surfaces influencing their possible 
evolution. In particular, when discontinuity surfaces are structured, in the 
sense that they are endowed with surface energy, surface interactions may 
occur. They are of standard and substructural nature. Standard surface 
stresses are in fact due to the deformation of the discontinuity surface while 
the other are a consequence of changes of substructure at the surface itself. 

Surface balances involving standard surface stresses have been obtained 
in [17] while substructural measures of interactions have been introduced 
in [22], [23] and the relevant balance equations derived there (see also [8]). 
However, there is no proof of their covariance. Such a proof is provided 
here (see Theorem 2) and is the main result of this paper. In obtaining it 
we enlarge the standard notion of observer: for us, in fact, such a notion 
involves not only the representation of the ambient space and the time scale, 
but also the representation of the manifold of substructural morphologies 
(see also [25]). 



2 Morphology of complex bodies 

We consider a body occupying a regular region £>o of the three-dimensional 
Euclidean space £ 3 (with affine translation space Vec) 1 . The current mor- 
phology of the body is described by two sufficiently smooth mappings: 

B 3X^x = x(X)ef 3 , B 3X^v = l>(X) eM. (1) 

1. x shows the current placement of a material element at X in Bq, is 
injective and orientation preserving, and B = x (So) is a lso regular as 
Bo- 



1 With the adjective 'regular' we refer to an open (bounded) subset of £ 3 co- 
inciding with the interior of its closure, with a surface-like boundary where the 
outward unit normal n is well-defined everywhere to within a finite number of cor- 
ners and edges. The treatment of infinite bodies requires only some minor technical 
adjustments in the results of the present paper. In any case, some remarks about 
them are presented throughout the paper. 
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2. i> is the order 'parameter map and assigns to each material element a 
coarse grained morphological descriptor u (X) of its substructure (or- 
der parameter), chosen as an element of a differentiable paracompact 
manifold M. (generally without boundary). 

Two natural tangent maps arise, namely Tx : TBq — > TB and Tv : 
TBo — * TM.. The pairs (x, Vx), (i/, Vi/) are the peculiar elements of Tx 
and Tf respectively. Since T£>o is a trivial bundle and a connection is natu- 
ral over it, we can separate in invariant way x from Vx which is commonly 
indicated with F. The condition that x be orientation preserving implies 
that, at each X G Bo, F G Horn (TxBo,T%rx.)B) has positive determinant. 
Since M. is a priori not trivial, the pair (i/, Vi/) cannot be separated in 
invariant way, unless there is a parallelism over M.. In principle, one may 
define in abstract way a parallelism over M. but, since the pair (i/, Vi/) en- 
ters constitutive issues, one should have a physically significant parallelism. 
In other words, when we act separating in invariant way Vi/ from v we 
should presume to have at least one physically significant parallelism over 
M. even when circumstances may allow us the use of invariance require- 
ments with respect to the choice of the connection. In fact, the presence of 
the pair (y, Vf) or of arbitrarily one element of it (namely v or Vi/) in the 
list of entries of the energy changes the representation of interactions. 

As remarked in the introduction, in dealing with the geometrical (mor- 
phological) description of complex bodies, we relax one of the axioms of the 
mechanics of simple materials. In our picture, in fact, the material element 
is not morphologically equivalent to a 'monad', a simple material particle 
in the sense of Cauchy (see [29], [30]), identified only by its place in space. 
Instead, we consider the body as a collection of subsystems of the same na- 
ture (the material substructure) and the order parameter at a given point 
represents the characteristic features of the morphology of the subsystem 
there. 

Sometimes the material substructure is a perfectly identifiable Lagrangian 
system as in the case of ncmatic liquid crystals [14], [20], [11] (in which the 
ncmatic molecules can be separated from the melt), sometimes it does not as 
in granular gases [6] and microcracked bodies [26]. In granular gases, e.g., a 
material element collects a family of sparse granules with peculiar velocities, 
so that the order parameter may be an element of a suitable Grassmanian 
of the tangent bundle of some finite-dimensional manifold, while for micro- 
cracked bodies each microcrack can be considered either as a sharp planar 
defect not interpenetrated by interatomic bonds or as an elliptic void, so it 
does not exist per se and has a volatile substance determined just by the 
surrounding matter. 

We leave undetermined the specific nature of u to cover a class as large 
as possible of special cases, following in this way the unifying point of view 
of Capriz [4]. Our primary strategy is to work with minimal requirements 
for M. and to add geometrical structure to it only when necessary. 
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In this way we avoid the path leading to a monotonous zoo of special 
models with slightly different formal aspects but with the same intrinsic 
nature. 

With respect to Hamiltonian elasticity of simple bodies (see [27], Chap. 
4) , here the main source of difficulties is the circumstance that M. is a non- 
trivial manifold; in particular it does not coincide with a linear space in 
most cases of prominent interest. 

The space C of pairs of maps (x, v), a product space of the type C x xC v , 
with x pertaining to C x and v to C v , has a non-trivial structure which 
depends on the geometrical properties of M.. Basically, we imagine that C x Q 
W 1,p (Bo, Vec) for some p > 1 and C v — PC 1 (Bo, M), i.e. we require that 
the order parameter map be at least continuous and picccwise continuously 
differentiable over Bo, while the placement map be an element of the Sobolcv 
space W ' p (Bo, Vec), even if we may basically require that x be continuous 
and piecewise continuously differentiable too. 

Specific examples showing the possible intricate nature of C, in particular 
of C u , can be discussed at length. 

We analyze below the case in which M. has a Riemannian structure. 
However, the developments in subsequent sections do not require strictly 
such a structure that will be called upon only when necessary. 

In any case, even when M. is Riemannian, we do not assign a priori any 
prevalent role to Levi-Civita connection. 



2.1 The case in which M. has a Riemannian structure: aspects of the 
nature of C 

We assume just in this section that M. is a Riemannian, connected, finite 
dimensional manifold, endowed with a metric gM an d a consequent Levi- 
Civita connection. 

Notice that we do not require here the embedding of M. in some linear 
space. When, in fact, we embed M. in a linear space by using, let say, 
Nash's isometric embedding to preserve at least he quadratic part of the 
kinetic energy (if it exists) pertaining to the substructure, the embedding is 
not unique. So the choice of it becomes matter of modeling. Moreover, the 
embedding is also not strictly necessary to build up the structures needed 
for the basic aspects of mechanics. For this reason, with the aim to eliminate 
overstructures as much as possible, we consider M. in its abstract setting, 
not embedded a priori in a linear space. 

We denote with (-,-) TM the scalar product over T V M. associated with 
gM- For any C^-curve [0, s*] 3 s i — > v (s) G Ai, we have a vector field 
v (s) = 7p (s), and for any other vector field A (s) over M. one defines its 
derivative along s \ — ► v(s) writing 2£- = V„( s )j4, and call A autoparallel 
along s i — * u (s) when ^j- = for any s £ [0, s*]. As usual, s i — ► v (s) is 
a geodesic when s i — > v (s) is autoparallel along it. 



6 Chiara de Fabritiis, Paolo Maria Mariano 

The Ricmannian structure assures the existence of a natural distance 
d_\4 : M. x M. — > M + over M.. In fact, by indicating with A any arbitrary 
piecewise C 1 curve A : [0, s*] — > A1 such that A (0) = v\ and A(s*) = i/ 2 , 
as usual we put 

d_M(i/i,i/ 2 ) = iiif {*(*)}, (2) 

varying A between V\ and 1/2, with I (A) the length 

(see, e.g., [1] Chapter 2, [33] Chapter 4 ). Locally, cLm is calculated over 
geodesies. 

However, cLm may be unbounded. In this case, we select a new metric 
d M : M x M — > M + over .M defined by 

dvvt 



J.M 



(4) 



1 + rfx 
Notice that the metric cLm 

(a) is equivalent to <1m (it leads to the same topology), 

(b) is bounded, 

(c) is of class C°°(M), 

(d) if cLm is complete on A4, then dj^i is complete too. 

A well know theorem due to Hopf and Rinow ([19] , p. 56) asserts that the 
following statements are equivalent: (i) dja is complete, (ii) closed bounded 
sets with respect to dj^t are compact, (iii) geodesies for dj^\ can be contin- 
ued as to be defined on the whole real line. Moreover, if M. is complete, 
the distance dj^i can be computed over geodesies not only locally, but also 
globally. 

Remark 1. When M. is compact, dj^ : A4 x M. — > M + is bounded, then 
it is not necessary to substitute it with dj^. Of course, the definition of dj^ 
is just one of the possible ways to bound the metric (i.e., the 'diameter' of 
the manifold is finite). 

By making use of dj^ we may define three natural distances over C v . 

1. The first metric distance, indicated with d^ l \ is defined by 

d^{v 1 ,v 2 )= J d M (u l7 u 2 )d(vol), (5) 

J Ba 

where (i) stands for 'integral'. 

2. Let {K n } be an exhaustion of Bo, i.e. a compact cover of So such that 

o 

K n C K n+ \ for any n. The second distance, indicated with d^ c ' , is 
defined by 

d (c) (£>i,£> 2 ) = y"^r max d M {v^u 2 ), (6) 

riGN 

where (c) stands for 'compact'. 
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3. More simply, a third natural distance, indicated with d^' , can be defined 

b y 

d (s) (£>i,£>2) = sup d M (1/1,1/3), (7) 

xeso 

where (s) stands for 'suprcmum'. 

Remark 2. There is a peculiar physical difference between the metric 
d^ and the other two. In fact, the exhaustion of Bo by compact sets, used 
to define d^ c \ implies that the values of the order parameter map over 
boundary points of £> may contribute slightly to the distance in C v , on 
the contrary of the other two metrics. In a certain sense, d( c > seems to be 
preferable when there are uncertainties in the physical meaning of boundary 
data about u, a problem (the one of boundary data) that appears very subtle 
for some material substructures like microcracks. 

Remark 3. Notice that d^ and d^ are complete over C° (Bo,M) if 
dj^t is complete. In particular, compactness of M. is sufficient for the validity 
of this statement. 

Remark 4. We realize that the space of continuous maps between Bo 
and M., namely C (Bq,A4), may be in general not complete with respect 
to the metric rfw. For example, let Bo coincide with the cube (—1, 2) in a 
frame OX1X2X3. We may analyze two different cases. 

(i) Case 1: M is not compact. Let us assume M = R, and for any x, y 6 K, 
put dM (x,y) = \x — y\. Let us construct also a sequence of v such that 
for (X 2 ,X 3 ) ranging in (-1,2) 

[ , if - 1 < Xi < 
u n = /„ (X) = \X? , if < Xx < 1 . (8) 

1 , if 1 < Xx < 2 



From J21, for ra < m G N, we get 



d (l) (fn, fm) = / d M (fn (X) , f m (X)) d (vol) < 

< f d M (f n (X)J m (X))d(vol) = 

J Bo 

1 1 



n +1 r7i+l 



0) 



so that the sequence /„ is Cauchy's. However, if we let n to infinity, the 
resulting limit is 

j (X) " \ 1 , if 1 < Xx < 2 ' (10) 

which is not continuous. The physical situation covered by this example 
is the one of a two-phase material in which the order parameter is the 
void volume fraction of one of the two species. The limiting distribu- 
tion l|10(l describes the complete separation of the two species with the 
occurrence of a discontinuity surface in the plane X\ = 1. 
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(ii) Case 2: A4 is compact. Let us assume Ai = S 1 C C and define with a 
slight abuse of notation 

[0 , if - 1< X x < 

u n = /„ (X) = \ exp (iX?) , »/ < Xi < 1 , (11) 

[ exp (i) , ifl<Xx<2 

for (X2, X3) ranging in (—1, 2) , with i the imaginary unit, and d,M (a, b) = 
oIm (Qj b) = \a — b\. For n < m, we get 

d (i) (/», / m ) = / dM (/„ (X) , f m (X))d (vol) = 

-'Bo 

= 9 / ^2 (1 - cos (Xf - X™))dXi < 

< 9 / (Xf - AT*) dXx = 9 (-L - -i-) , (12) 
Jo \n + l m + lj 

so that the sequence /„ is Cauchy's but, as n — > +00, its limit /(X), 
with 

; W ~ \ exp (») , i/ 1 < Xi < 2 ' (W) 

is not continuous. The physical situation covered by such an example 
could be (rather ideally) the one of a porous medium with connected 
voids full of normal supernuid Helium (i.e., 4 He; see [4] pag. 7) which, to 
the limit lfH3)l . percolates on one side with respect to the plane X\ = 1. 

Remark 5. The metric (JHJ may furnish unreasonable results for un- 
bounded bodies. In fact, let us consider a Cartesian frame OX1X2X3, and 
a domain B = D q x [0, +00) where D q is the square (—1, 1) in the plane 
X3 = 0. Imagine that the closure of the domain B is occupied by a semi- 
infinite beam made of a two-phase material so that the order parameter is a 
scalar ranging in [0, 1] and representing the volume fraction of one phase. We 
consider a distribution of the two phases homogeneous with respect to X3 
so that in a given configuration the distribution of the phases is expressed 
by the order parameter mapping B 3 (Xi, X2,X^) — ► vx = i>x (^1^2) £ 
M. = [0,1]. We consider another distribution of the two phases, namely 
v% = 1 — vx and select oIm {x, y) = d,M (x, y) = \x — y\ so that the integral 
Jn dM i^x, V"i) d (area) = J D \2vx — 1| d (area) is a certain constant posi- 
tive K (when v\ is not a constant equal to 5), then the distance S- 1 ' (vx, ^2) 
reduces to the integral of K over [0, +00) which is +00. In other words, 
the two configurations obtained through a rigid rotation of the beam are 
at infinite distance if one uses the integral distance (JjJJ. The same result 
does not occur with the distances © and J7J) and the calculation reduces 
basically to the evaluation of the maximum of the function |2z>i — 1| varying 
(Xx,X 2 )vaD q . 
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To control the behavior of the derivatives of v , basically of Vz> with 
Vf = VJ> (X) G Hom(TxBo,T^A4), we find problematic to act directly 
on Hom(TBo,TA4). In fact, when we select i>\ and £2, at a given X, we 
get in general V\ ^ i/ 2 so that we get Vi'i £ Hom{TyJ3 a ,T Vl M) and 
V1/2 6 Hom(TxBo,T V2 A4). Though g^vi induces natural metric structures 
over T* M. and tensor product bundles of TM. and T*A4, it does not give a 
natural way of comparing V1/1 and V1/2 because Vf is in certain sense like 
the two-point tensor F. To compare them it is then necessary to transport 
\7v>2 over Ui through a connection or Vi>i over 1/2, i.e. to transfer elements of 
Tl, 2 M over T Vl M.. Since we do not assign any prevalent role to Levi-Civita 
connection, even in the case in which M. is complete and we may connect 
i>\ and 1^2 with a geodetic, we face the circumstance that the transport is in 
general not isometric, so that the comparison (defined in some way) of the 
two derivatives at V\ may lead to a different result of the same comparison 
at V2- Moreover, even in the case in which one would choose to use both the 
comparison at u 1 and the one at 1^2, say, using half of their sum, one could 
face the basic difficulty that the transport could be unbounded (due, e.g., 
to topological features of M) or the bound could not be uniform over the 
curve connecting V\ with U2- In this case, the possible choice to consider 
admissible only values over A4 that can be connected by curves assuring a 
uniformly bounded transport could reduce too much the generality of M. . 

By considering previous remarks, we follow a procedure leading to the 
possibility to compare Vvi and Vi^ in Horn {T%Bq,T^Bq). 

To this aim, we consider first the adjoint of Vf^, indicated with \Ji>* and 
such that Vu* = VJ>* (X) e Horn (T* ± M, T^B ), then we define 

J a< (Vi/ 1 ,Vi/2) = ||Vi^Vi/i-Vi/5Vi/2|| (14) 

where here ||-|| is the usual norm in M 3 ® M 3 , i.e. the usual norm of 3 x 3 
matrices such that, for any A e M 3 <g> M 3 , one has || A|| = ^Jtr (A T A). 

Remark 6. We express Vi>*Vu as {Vu*Vv) AB — Vi>*^ (g w ) « Vi>^ 
in components. In other words, the second order tensor V£'*Vt' is the 
pull-back in Bo of the metric gx in A4. As a consequence, the distance 
d-M (Vi'i, V1/2) compares values of the metric gM at two different points 
(namely V\ and 1/2) of M.. In fact, the scalar product Vi/*Vi'i • (dX. <8> dX) 
is the 'length' of dv at v\. 

With d,M as above, we define the counterparts of the distances JSJ , I© 
and J7|. 

4. The first distance, indicated with S l \ is defined by 

J w (Vi>i,V£>2)= / (f M (V^i,Vi/ 2 )d(wO. 



(15) 



5. Let {K n } be an exhaustion of Bo, i.e. a compact cover of Bq such that 

o 

K n C K n+ i for any n. The second distance, indicated with dS c \ is 
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defined by 

d^ (Vi>i,Vi>2) - Y, ^S dM ( Vvi > Vv *) ' ( 16 ) 

n£N 

6. More simply, a third natural distance, indicated with dS s > , can be defined 

by 

d (s) (Vi>i, Vi> 2 ) = sup d M (Vi/i, Vi/ 2 ) , (17) 

xeBo 

when sup XeSo d>i remains bounded over Bo- 

Remark 7. We do not introduce any 'normalization' of the distance like 
the one in Q. In fact, (|15|) - (|17fl are calculated over bundles whose fibers 
are linear spaces so that dj^4 (•, •) displays possible properties of homogeneity 
(natural over fiber spaces) while a normalization like Q would not. 

Appropriate topologies on C v may then be induced by the distances 

~d (i) (£>i,£> 2 ) = d (l) (i>i,i>2) + d w (Ve>i, Vi> 2 ) , (18) 

<? c > C£>!,£> 2 ) = d^ (£>!,£> 2 ) + J^ (Vf>i, V£ 2 ) , (19) 

rf (s) (£>i,£> 2 ) - d w (i>i,i>2) + d (8) (Vi>i, Vi> 2 ) . (20) 

Remark 8. A program based on the systematic analysis of the proper- 
ties of C v associated with each metric above could deserve to be followed. 
In fact one could analyze systematically the influence of the topological 
properties of M. (such as compactness) on the structure of C v . Here, such 
a program goes far from the main purposes of the present paper and we 
do not follow it. However, the remarks above underline the nature of subtle 
geometrical difficulties that one may encounter in going on and constitute a 
reasonable starting point for further investigations. In any case, one should 
take into account that when the manifold M. is per se a linear space or for 
reasons of modeling is embedded isometrically in a linear space, the results 
in [2] and [18] apply directly to characterize the topological properties of 
the space C. 

2.2 Something more about kinematics 

Motions are sufficiently smooth curves over C. For a given interval of time 
[0,t], we then have mappings [0, i\ 3 t i — > (x t ,/> t ) G C and indicate with 
x = x (X, t) and v = v (X, t) the current place at t of a material element 
resting at X when t = and the current value of the order parameter. 

With x and i> we denote rates given by 4f (X, t) and 4£ (X, t) respec- 
tively, with v G Tj>(x.t).M. They have counterparts v and v in the current 
place B given by 

6x[0,(]3(x,()^v = v(x ! t)eT I 6 (21) 
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and 

B x [0, i\ 9 (x, t) h£+ v = v (x, t) S r i/(x(Xl t),t)-M, (22) 

obtained through the mapping X i — > x = x (X) at each t. We have x = v 
but v = i> + (gradv) v and we may write also v = v + (Vf) F _1 v = i> 
(Vis) X, where X is the material velocity — F _1 v associated with the inverse 
mapping X = x (x, t). 

For the acceleration of the order parameter we have 

i) a = d t v a + r^^ (23) 

with 

r^ = \sm (^ (%M) 0S + d v „ (g M ) 7 * - d v s (g M ) 0J ) (24) 

Christoffel symbols relevant for M.. The expression (I23|l enters the repre- 
sentation of possible inertial terms pertaining to the substructure 2 . 

Let g be the metric in the ambient space (in general, for curved frames 
we have x i — ► g (x), i.e. g depends on the place), by indicating with F T 
the transpose of F, the mapping 

6o3X^ (F T F) = C (X) e Sym+ (T X S , T^B Q ) (25) 



is the pull-back of g through the deformation x and in coordinates we have 
C AB ={F T ) A i g ij Fi. 

If Bo is endowed with a not flat 'material' metric 7 (i.e. Bq 9 X h^% 
7 = 7 (X) G Sym + (Vec, Vec)), the difference (C — 7) (X) is twice the non- 
linear deformation tensor E (X) which measures relative changes of lengths 
by using Bq as paragon setting. In an alternative point of view, we may 
consider B as paragon setting pushing forward 7 and comparing lengths 
there as explained in all treatises on non-linear elasticity of simple bodies 
in chapters dealing with measures of deformation. 

In the case of complex bodies the matter may be more intricate and 
a general (in certain sense abstract) treatment of measures of deformation 
seems to be absent (the only exception being the basic remarks in [5]). 
The key point is the specific nature of v. In fact, when u represents, e.g., a 
microdisplacement, an independent rotation or an independent deformation, 
its gradient enters the measures of deformation together with v itself. On the 
contrary, when v describes a property not related strictly with changes of 
lengths (say in the case in which v represents the volume fraction of a phase 
in a two-phase material or the spontaneous polarization in ferroelectrics 
etc.), standard C and E, or their spatial counterparts, are sufficient to 
measure the macroscopic deformation. 

In general, we could imagine to have a map of the type 

(F, g, i/,Vi/) ^ G (F, g, i/.Vi/) 6 Sym+ (T X B , T x 2? ) (26) 



Basic remarks about the influence of f on inertia are discussed in [7]. 
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with G a metric in Bq involving the pull-back of g, and to define a defor- 
mation tensor E (X) as | (G — 7) (X), erasing v and/or Vf each time in 
which circumstances suggest such a cancellation. 

Of course, the information furnished by 1)26(1 is rather volatile unless we 
put it in a context specifying its nature for some particular substructure. 

Remark 9. Cosserat [10] and micromorphic ([28]) materials are well 
known classes of complex bodies in which the order parameter (or its gra- 
dient) affects the representation of the measures of deformation (see also 
essential remarks in [4]). In the former class each material element is pic- 
tured as a 'small' rigid body that may undergo rotations independent of 
the surrounding material elements. Of course, neighboring relative rota- 
tions may alter lengths, in a sense described e.g. in [10] or [4]. In the latter 
class one imagines that each material element may suffer independent de- 
formations (it is like a ball of rubber). Measures of relative deformations 
may be then introduced. Less popular is the case of microcracked bodies 
and we give some details about it in order to construct an explicit example 
of the possible influence of the order parameter and its gradient on the mea- 
sures of deformation. In fact, when microcracks are smeared throughout a 
body, the material element is pictured as a 'patch' of matter endowed with 
a population of microcracks that can be considered either as planar sharp 
defects not interpenetrated by interatomic bonds or as elliptic voids with 
non-null volume and one dimension very small with respect to the others. 
If we consider frozen the microcracks in a given material element placed at 
X in £>o, a standard deformation x puts it (or better its centre of mass) 
in the place x = x(X). Now, if we allow the microcracks to deform (say 
without growing irreversibly for the sake of simplicity), the centre of mass 
of each material element undergoes in principle a shift 3 toward a new place 
x'. If we indicate with B' the minimal regular region containing the col- 
lection of x', each one corresponding to each x, we may imagine (this is 
a basic point of modeling) to obtain B' from B by means of a sufficiently 
smooth mapping f such that x' = (f o x) (X) and f (£>) = £>'. By denoting 
with grad the gradient with respect to x, as before, by chain rule we get 
V (f o X) (X) = ((gradf) Vx) (X) = F( m )F, where F( m ) = (gradf) (x) is the 
gradient of deformation from B to £', i.e., F^ = Horn (T X B, T*B'). If we 
indicate with d (X) = (fox) (X) — x (X) the displacement from x to x', 
defined as a field over Bq, since Vd (X) = (gradd a )F, with d a = do x _1 , 
we get the additive decomposition V (f o x) (X) = F+Vd (X). If we write 
F to t — F+Vd (X), F to t is the gradient of deformation from Bo to B' . So that 
the right Cauchy-Green tensor Ctot (X) = Ff ot F to t involves Vd. Moreover, 
by comparing the additive decomposition of F to t with the multiplicative 
one, namely F^F, we realize that F^ = I+Vd (X) F" 1 = I + gradd a . 
The direct (perhaps naive) description of the kinematics of microcracked 



3 Here the circumstance that there is a population of microcracks in the material 
element and that each microcrack is not a spheric void is crucial. 
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bodies just sketched here follows [21], [26]; however, it can be obtained by 
using the procedure involving the limit of bodies described in [12], [13]. 

2.3 Observers 

Three geometrical environments are necessary to describe a motion and 
then the configuration at a prescribed instant of a complex body: the time 
interval [0,t\ (that we could extend in principle to the whole real line), the 
ambient space £ 3 (of course £ 2 or £ l if we deal with two- or one-dimensional 
bodies), and the manifold M.. 

In classical mechanics an observer is a representation of the measure 
of time and the ambient space so that changes of observers are changes of 
representation. 

Here, our point of view is that the notion of observer should involve 
all the descriptors of the morphology (i.e. the geometry) of the material 
element. Then, for us, an observer O is a representation of a triplet of 
geometrical settings: (i) the interval of time, (ii) the ambient space £ 3 and 
(iii) the manifold M. of substructural states [25]. Such a definition is a 
natural extension of the standard point of view including just (i) and (ii): it 
takes into account the presence of the manifold of substructural states M. 
as natural geometric ingredient in the representation of the body. 

For the sake of simplicity we shall consider observers coinciding about 
the measure of time. 

General changes of observers, as used below, will then involve automor- 
phism of £ 3 and the action of arbitrary Lie groups over M. In this sense 
the notion of change of observers involves a pair of one-parameter families 
of transformations indicated below. 

Al. R+ 3 s 2 i — > f 2 2 £ Aut (£ 3 ), with f 2 the identity 4 . We put f 2 ' (X) = v, 
where the prime denotes differentiation with respect to S2- 

A2. A Lie group G, with Lie algebra g, acts over M. If £ € 0, its action 
over v &A4 is indicated with t; M (y). By indicating with v g the value 
of i> after the action of g £ G, if we consider a one-parameter smooth 
curve K+ 3 s 3 i — ► g S3 6 G over G, such that £ = -Jj 2 - | S3= o, and 
its corresponding curve S3 1 — ► v 9s over M., starting from a given u, 
we have £\ M (u) — -r~v ga | S3 =o- When G coincides with SO (3), then 
t; M (v) = A (v) q, with A (is) defined above and qA € so (3). 

We specify for future use in special cases the relations characterizing 
'rigid' changes of observers, i.e. the ones characterized by the contemporary 
action of SO (3) on M. and £ 3 . The latter action is induced by standard 
isometries of point space. 

Let [0,t] 9 t 1 — > Q(t) e SO (3), with Q (0) =Id SO (3), b e a smooth 
curve on SO (3). Let also O be an assigned representation of £ 3 at t = 0. 

4 Aut (£ 3 ) is the group of automorphysms of £ 3 . 
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At any t ^0 we may associate an observer 0' t obtained from O by means 
of Q (t). The transformation from O to 0' t is isometric, then a point x seen 
by O is mapped in a point x' = w (t) + Q (t) (x — xo), where t i — > w (£) 
is an arbitrary point valued function smooth in time and Xo a fixed point 
chosen at will in space. If we calculate the rate of x' and pull-back it in O 
writing x* for Q T x', we get the standard reaction 

x* = x + c (t) + qA (x - x ) , (27) 

characterizing a classical change of observer, where c = Q T w is the transla- 
tional velocity and q the rotational one. If the arbitrary element qA of the 
Lie algebra so (3) acts also on A4, it induces a rate (that we call 'rigid', with 
the subscript R, just for reminding the circumstance that it is related to a 
spatial rigid body motion) given by „4q, where A(y) £ Horn (Vec, T U A4). 
Classical changes of observer in a multifield setting are then characterized 
by the relations (|27[l and 

v* = i> + Aq (28) 

when they agree about the measure of time. Basically, if v q denotes the value 
of v after the right action of SO (3) over Ai, we have A = -j 3 - | q= o , where q 
is a vector connected with Q £ SO (3) by the formula Q = exp (— eq).with 
e Ricci's permutation index. For example, let M. be coincident with S 2 as in 
the case of magnetostrictive materials). For any t £ 5 2 , we get x q = Qt, 
then A = — rA. 



2.4 Relabeling 

Another transformation playing a role in the developments presented in 
what follows is the relabeling of material elements in the reference place 
Bo. With the term 'relabeling' we indicate a generic transformation of £>o 
induced by a C 1 (Bo) point valued mapping f 1 such that (i) it is isocoric 
(i.e. it preserves volume) and (ii) establishes also a diffcomorphism between 
Bo and f 1 (Bo). From a physical point of view, the important part of the 
action of f 1 relies in a sort of 'permutation' of possible defects or, better, 
inhomogeneities in Bo- In particular, we will consider a one-parameter family 
of such isocoric diffeomorphisms defined below: 

A3. R + 3 si 1 — ► f s \ £ SDiff(Bo), with fg 1 the identity; i.e. at each si 
we get X 1 — ► f (X), with Divfj' (X) = 0, where the prime denotes 
differentiation with respect to the parameter s\. We put fg' (X) = w. 



3 Langrangian 3+1 forms and balance equations 

The multifield theoretical analogue of the theory of elasticity for complex 
bodies (i.e., bodies with material substructure) relics on rather articulated 
fiber bundles. 
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We start by considering a fiber bundle it : y -^Bq x [0, t\ such that 
7r _1 (X,t) — £ 3 x M. A generic section r\ e r (y) of y is then i] (X,t) = 
(X, i, x, v). For sufficiently smooth sections, the first jet bundle J 1 ^ over y 
is such that J x y 3 j 1 (r?) (X,£) = (X, t, x, x, F, u, 0,\7v>). 

Up to this point the discussion has been purely geometric. No issues 
related with the constitutive nature of the body and with the interactions 
arising inside it have been discussed. They enter into play here: we assume 
that the general body under examination is made of (non-linear) elastic 
material. In this case we may associate with it the canonical Lagrangian 
3 + 1 form 

L: J 1 ? -► A 3 * 1 (Bo x [0,i\). (29) 

The definition of the elements of the space of 3 + 1 forms A 3+1 (Bq x [0, t\) 
would require some care. In fact, Bo X [0, t\ is a manifold with boundary 
coinciding only with {0} x Bq U {£} x Bq (because Bo is open and coincides 
with the interior of its closure) , while the definition of odd forms is imme- 
diate on manifolds without boundary. However, one is commonly interested 
in evaluating the variation of the total Lagrangian 

Z(B )= f £d 3 XAdt; (30) 

JB x[0,i\ 

so that, in defining L (j 1 (77)), possible problems related with boundary 
points do not play any role. 

We assume that L admits a density C such that 



L (j 1 (77) (X,t)) = C (X, x, x, F, u, u,Vu) d 3 XAdt (31) 



with 



1 2 

C (X, x, x, F, v, u,Vv) = -p |x| +p oX (v, v)-p e (X, F, v,X7v)-p Q w (x, v) , 

(32) 
where p is the referential mass density (conserved during the motion), \ 
the substructural kinetic co-energy (needed when physical circumstances 
prescribe substructural inertia), e the elastic energy density and w the den- 
sity of the potential of external actions, all per unit mass. The presence 
of u in the list of entries of x an d e is due to the non-trivial structure of 
Ai. In general the elements of TA4 cannot be separated invariantly unless 
a parallelism can be found over M.. From classical non- linear field theories 
we know that x disappears from the list of entries of e due to reasons of 
invariance, but v does not (see [4]). 

We assume that C (•) be sufficiently smooth so that we may find at least 
one section rj satisfying Euler-Lagrange equations 

&Z = d x C - DivdpC, (33) 



d C = d u C - Divd v „£. (34) 
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With respect to the families of transformations characterizing relabeling 
and changes of observers, for the sake of brevity, we indicate with f , f 2 and 
v g the values f s \ (X), f 2 2 (x), v B , a (X). 

Definition 1 A Lagrangian density C is invariant with respect to the action 
off^fl and G if 

£(X 1 x,4,F,i/,i>,Vi/) = 



= C (f\ f 2 , {gradi 2 ) x, (gradf 2 ) F (Vf 1 ) \u g ,u g , (Vv g ) (Vf 1 )" 

'(35) 
for any g 6 G and Si,32 € K + - 

Let Q and 3 be scalar and vector densities given respectively by 

Q = diC- (v - Fw) + di,C- {i M (u) - (Vi/) w) , (36) 



$ = £w+ {d F C) T (v - Fw) + (dv.Cf {Z M W) - (Vi/) w) . (37) 



Theorem 1 [9] . If the Lagrangian density C is invariant under f^ , f 2 2 
and G, for any g £ G and Si, S% S K + , then 

Q + Div$ = 0. (38) 



The detailed proof is contained in [9]. Here we remind only that it is 
based on a direct explicit calculation of the terms in l|38|l on the basis of 
(|36|l and l|37|) , and on the exploitation of the relations 



_d_„, d 

is\ 
d 



-C | Sl= 0.S2=0,S 3 =0 — 0, C | Sl= 0,S2=0,S3=0 — 0, 

ds\ ds2 



-C | Sl= o.s 2 =o,s 3 =o —0, (39) 

ds 3 

which are consequences of the request of invariance for C. 

Corollaries [9] . The following statements hold true for any Lagrangian 
density of the type 1-il'l) in the sense of Definition 1: 

1. If f 2 2 alone acts on C leaving v arbitrary, from (|38|l we get Cauchy's 
balance of momentum 

p x =p Q b+DivP, (40) 

where P = -d F £ e Horn (T£S , T*B) is the first Piola-Kirchhoff stress 
and b = <9 X £ € T*B the vector of body forces. 
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2. If G arbitrary acts alone on C, from 1(38(1 and the arbitrariness of the 
element £ of the Lie algebra of G (chosen to define ^ M (u)) we get 
Capriz's balance of substructural interactions 



Po [9uX - 9uX) = -* + Po/3 m + DivS, (41) 

in covariant way, where (3 = —p^d u w G T*M. represents bulk non- 
inertial external interactions acting on the substructure, S = —d^j v L £ 
Horn (T^Bq, T*M) takes into account contact substructural interactions 
between neighboring material elements (and is called microstress) and 
z = —p d v e £ T*J\A indicates self- interactions of the substructure in 
each material element (and is called self-force). 

3. Let G = SO (3) and, for any element qA of its Lie algebra, f? be such 
that v = qA(x-xo) with Xo a fixed point in space. In other words, we 
require that the same copy of SO (3) acts both on the ambient space 
and on M. If we require the invariance of e with respect to SO (3), we 
obtain 

skw (d F eF T ) = e (a t d v e + (V^ T )' d w „e\ , (42) 

where e is Ricci's alternator and skw (■) extracts the skew-symmetric 
part of its argument. Previous statement render more perspicuous Re- 
mark 3 of [9]. 

4. If f * alone acts on C, with w arbitrary, from 1)38(1 we get 

(F T di£ + ^v T d C) - Div (f- Qp |x| 2 + PaX (v, *) J i) - M = 

(43) 
where P = p el - F T P - Vv T S e Aut {T£_B ), with I the second order 
unit tensor, is the modified Eshelby tensor for continua with substruc- 
ture proposed in the general setting in [22], [23]. In particular, (|43|) co- 
incides with the balance of configurational forces for a continuum with 
substructure in absence of dissipative internal forces driving defects. The 
balance ((43(1 . in fact, is only a consequence of the invariance with respect 
to relabeling. 

5. Let G = SO (3) and, for any element qA of its Lie algebra, f^ is such 
that w =q A (X — Xo) with Xo a fixed point in Bq. If the material is 
homogeneous, and only the special choices of f£, and G just defined act 
on £, P is symmetric. 

Remark 9. It may be asked whether or not an integral version of 1(41(1 
can be postulated as integral balance principle of substructural interactions 
and then used as a first principle. In general the answer is negative unless M. 
is a linear space or is embedded in a linear space. The reason relies upon the 
circumstance that the eventual integrands /3, z and <Sn would take values in 
T* M. which is a non-linear space unless the above mentioned situations of 
linearity for M. occur, so the integrals are not defined. In other words, each 
time we use an integral version of 1(4111 we presume implicitly an embedding 
in a linear space; on the contrary it does not make sense. 
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4 Discontinuity surfaces 

In common special cases, solutions of <|33ll . (13411 or, with other notations, of 
(|40[1 . (|41|l are not smooth and may display discontinuities concentrated on 
submanifolds of codimension 1. Moreover, experiments display domain for- 
mation and branching of microstructures of various nature (see e.g. cases of 
nematic order in liquid crystals, polarization in ferroelectrics, magnetization 
in micromagnetics, superconducting domains in superconductors etc.). 

The presence of the gradient of the morphological descriptor in the list 
of entries of the energy accounts for the presence of interfaces of domain 
walls in a smeared sense. However, one may ask what happens when addi- 
tional discontinuity surfaces occur and there is interaction between smeared 
and sharp interfaces. As an example, the us consider the solidification of a 
two-phase flow [24]. Two types of interfaces occur in this case. Interfaces 
between the two phases of the fluid and the interface between solid and fluid 
parts. The two-phase fluid is a complex body, the material element is a sort 
of box containing the two phases in different portions, so that a natural 
morphological descriptor is the volume or mass fraction u of one phase with 
respect to the other. The presence of the gradient of u) in the list of entries of 
the energy allows us to account for diffuse interfaces between the two fluid 
phases. These diffuse interfaces interact with the interface solid-fluid. The 
latter may be described either either as a diffuse interface too or a sharp 
interface. The latter interpretation is a special case of the theory discussed 
below. 

Really, even in single phase complex materials the energy may depend 
on the gradient of the order parameter to account for the inhomogeneous 
behavior of substructures. For example, in the case of micromorphic mate- 
rials, each material element is a cell that may undergo micro-deformation 
independently of the neighboring cells, in addition to the participation to 
the overall macroscopic deformation [28]. The morphological descriptor is 
then a second-order symmetric tensor representing such an additional inde- 
pendent deformation. Such a situation may be representative of the behavior 
of elastomers, for example. Adjacent material elements may then undergo 
in principle different micro-deformations so that, in going from one element 
to the other, the energetic landscape changes in accord to the gradient of 
the micro-deformation and weakly non-local interaction effects of gradient 
type between neighboring material elements may be accounted for. Even 
in this case, additional discontinuity surfaces may occur as defects or shock 
and acceleration waves. Energetic effects associated with the gradient of the 
micro-deformation may influence them. 

Analogous situation occurs also in quasicrystals where collective atomic 
modes (phason degrees of freedom) occur within each crystalline cell (the 
material element) and are represented by a stretchable vector. The energy 
depends both on the standard strain and the gradient of such a vector [31]. 
Even in this case, discontinuity surfaces such as racks or dislocations may 
intervene. Their behavior is influenced by the energetic effects associated 
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with the gradient of the order parameter (here the vector mentioned above) 
but they are not strictly modeled through it. 

The physical examples above allow us to clarify the point of view followed 
below. We treat, in fact, situations in which smeared interfaces or gradient 
effects due to inhomogeneous behavior of material substuctures interact with 
sharp interfaces or, more generally, discontinuity surfaces. 

In particular, we focus our attention on a single discontinuity surface £ 
defined by 

£ = {XecJB, /(X)=0}, (44) 

with / a smooth function with non-singular gradient. It is oriented by the 
normal vector field £ 3 X i — > m (X) = V/ (X) / |V/ (X)| and we use the 
notation 77 for I m ® m. 

For any field e (•) continuous and piecewise continuously differentiable 
on £, we indicate with V^e its surface gradient at X, with e — e (X), given 
by V^e = Veil. The opposite of the surface gradient of m, namely — V^m 
is indicated with L and is the curvature tensor. 

Let X i — > a = a (X) be a generic field taking values in a linear space 
and suffering bounded discontinuities across £. For e > we indicate with 
a the limits lim E _*o a (X ± em) which are the outer (a + ) and inner (a~) 
traces of a at £. Then we denote with [a] = a + —a~ the jump of a across £ 
and with 2 (a) = a + + a~ its average there, so that for any pair of fields a\ 
and ai with the same properties of a we have [0102] = [ai] (02) + (ai) [02] 
if the product a\ai is defined in some way and is distributive. 

[F]U = implies that £ is coherent. In other words the two pieces of 
the body separated by the surface do not suffer relative shear. 

Moreover, if we attribute any 'virtual' motion to £ by means of a vector 
field £ 3 X 1 — > u (X) <E R 3 with normal component U = u • m and assume 
that the velocity x may suffer bounded jumps across £, we get the condition 
[x] = — U [F] m, whose proof is textbook affairs. 

As essential point we assume also that the order parameter map D is 
continuous across £. A special case in which such an assumption plays a 
prominent role is Landau's theory of phase transitions. 



Jf..l The unstructured case 

We treat first the case in which £ is unstructured, i.e. it is free of own 
surface energy. We derive the balance equations across £ in a way that 
prove their covariance which has been not discussed so far. The procedure 
we adopt relies on the exploitation of an integral version of the pointwise 
relation P%|) . 

We call part any regular subset b of £>o with non vanishing volume 
measure. We consider an arbitrary part bs crossing £ in a way in which 
the intersection of its boundary dbs with £ be a curve and write for it the 
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integral counterpart of 1)38(1 . namely 



^ f Qd (vol) + f $-nd 

at Jbj, Jabs 



(area) = 0. (45) 



If we take fixed the part bs with respect to the virtual motion of S, by 
the use of the transport and Gauss theorems (see, e.g., [27], [35]), as bs — > 
bs H U we find 



d_ 
db 



f Qd(vol)->- f [Q]Ud(vol) 

Jbs Jb s nE 



(46) 



/ £ • nd (area) -» / [$] • md (area) , (47) 

so that the arbitrariness of b.j; implies the pointwise balance 

-[G]tf+[3]-m = 0. (48) 

Proposition 1 // the transformations Al, A2, A3 are smooth throughout 
the body, the validity of \4$j) and the invariance of C imply covariant point- 
wise interfacial balances across S as in the list below. 

1. The action ofi^ 2 alone implies the interfacial balance of standard interactions 

[P] m = -p [x] U. (49) 

2. The action of G alone implies the interfacial balance of substructural 
interactions 

[S] m = -p [dox] U. (50) 

3. The action of f£ alone implies the interfacial configurational balance 
along the normal m in absence of dissipative forces driving S, namely 

m [P] m = Po U [(Vi/)* dox] -™+\p [x (", &)] ~ \p U 2 [|Fm| 2 ] . (51) 

The proof is a consequence of direct calculation of l|48|) . It is contained 
in the one of Theorem 2 below. 

Balance equations at discontinuity unstructured surfaces in standard 
elasticity of simple bodies are derived by means of a direct evaluation of the 
variation of the total Lagrangian in [15]. 



GEOMETRY OF INTERACTIONS IN COMPLEX BODIES 21 

4-2 The structured case 

We consider here S endowed with a surface energy density 4> associated 
with surface tensions of standard and substructural nature. 
Let F and N be defined by 

F = (F)77, N=(Vz/}77. (52) 

They are surface deformation gradient and surface gradient of the order 
parameter at X. There exist then two mappings F and N such that 

£ 3 Xh^F = F (X) e Horn (T*£, T X B) , (53) 



E 3 Xh^N =N (X) e Horn (T*£, T V M) . (54) 

Elementary algebra provides us 

(F) = F+((F)m)<g>m, (55) 

(Vi/)=N+((Vi/)m)<g>m. (56) 

The surface energy density is then defined by 

(m,F, v,N) i-^ </> = </> (m,F, i/,N) (57) 

and we assume that <j> be sufficiently smooth, ft is worth noting that the 
presence of the normal m in the list of entries of cf> points out that we are 
considering anisotropic surfaces. 

We require the invariance of (j> with respect (i) to general changes of ob- 
servers and (ii) to relabeling of S. As discussed above, changes of observers 
are characterized by the action of the group of automorphisms of £ 3 and of a 
generic Lie group over M. . However, definition A3 of f, 1 needs to be modified 
in order to describe the relabeling of S in addition to the overall relabeling 
of Bq. Instead of Si i — ► f* we should consider time-parametrized families 
si i — ► fg of elements of SDiff (Bo) characterized by the properties listed 
below. 

Relabeling of B§ including E. 

1. The map s\ i — ► f, satisfies Al. Moreover, we require that the field 
Bo 3 X i — > w = w (X) = fg (X) is at least of class C 1 (Bo), then also 
across and along S. 

2. Each f^ preserves the elements of area of E. Namely, if dA is the element 

of area of S in Bo, dA = f}* o dA, where the asterisk indicates push 
forward. 

3. (Vw)m = 0. 

4. V £W m = 0, with w m = w • m. 
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Two lemmas are useful for subsequent calculations. 

Lemma 1. For any isocoric vector field w (•) of class C 1 (Bo) , 

H ■ V r w = ((Vw) m) • m, (58) 

with w = w (X). 

The result follows from direct calculation. 

A second order tensor field S 3 X i — ► A = A (X) e Horn (R 3 ,R 3 ) is 
called superficial if Am = at each X. 

Lemma 2. For any second-order superficial tensor field A on S, one 
gets 

m • DivzA = A ■ L. (59) 



Definition 2. A surface energy density <j> is invariant with respect to 
the action of f* , ff 2 and G if 

4> (m,F, i/,N) = (Vf 1T m, (grad s f 2 ) ¥ (Vf 1 ) ~* , i/ fl ,N fl (Vf *) "*) , 

(60) 
for any g E G and Si,«2 € K- + , where N 3 = (Vi/ g ) II and we have used 
notations common to Definition 1. 

Let X be a sufficiently smooth vector density defined over S by 

X = -077w+(cW) T (v - (F) w) + (9n^) T (^ (i/) - (Vi/) w)-(9 m m) w. 

(61) 
It is the counterpart of 3" for a surface. 

Theorem 2. Lef £ be a structured surface with surface energy <j>. Let 
us assume 



d 



/ Qd (vol) + # • nd (area) + / X • nd (length) = (62) 



/or anj/ part bs of Bo crossing S. If C and (f> are invariant with respect to 
f} , ff 2 and G, covariant pointwise balances across U follow as in the list 
below. 

1. The action of i 2 alone implies the interfacial balance of standard inter- 
actions 

[P] m + Div s T = -p [x] U, (63) 

where T = —dw<t> G Horn (Tx.S, T£B) is the surface Piola-Kirchhoff stress. 
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2. The action of G alone implies the interfacial balance of substructural 
interactions 

[S] m + Div s S -$ = -/>o [do X ] U, (64) 

where §= —dn<fi G Hom(T^S,T*A4) is the surface microstress and 
3 = d u 4> £ T*M the surface self-force. 

3. The action of f, 1 alone implies the interfacial configurational balance 
along the normal m in absence of dissipative forces driving S, namely 

m- [P] m + C ta „ • L + Div s c = 

= Po U [(Vi/)* dox] ■ m + p Q {x (i/, 0)] - \p Q XJ 2 [|Fm| 2 ] , (65) 

where 

Ctan = (t>n- F T T - N T S (66) 

is a generalized version of the surface Eshelby stress and 

c = -d m <f) - T T (F) m - § T (Vi>) m (67) 

is a surface shear. 

By following different procedures, equation (|63[1 has been derived in [17] 
while equations l|64|l and l|65l) in [22], [23] (the version of l|65|) for simple 
bodies is clearly obtained in [16]). Here, with Theorem 2 we prove their 
covariance: this is the main novelty of the theorem itself. 

4-3 Proof of Theorem 2 

Step 1. 

In accord to Definition 2, since (f> is invariant, we have 



— 4> ( Vf 1T m, (grad s f 2 ) F (Vf *) , i/ ff) N 9 (Vf *) j | Sl=0 , S2 =o, S3 =o = 0, 

(68) 



■£-]> (vf 1T m, (gra^f 2 ) F (Vf *) * , i/ fl ,N, (Vf *) ^ \ Sl=0 ,s 2 =o,s 3 =o - 0, 

(69) 



^-J> f Vf 1T m, (gra^f 2 ) F (Vf 1 ) * , i/ g ,N s (vf 1 ) ^ | Sl=0 , S2 =o, S3 =o = 0, 

(70) 
which correspond respectively to 

F T T • V s w + N T § • Vsw + d m cj> ■ (Vw) m = 0, (71) 
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T-V s v = 0, (72) 

3-U(")+S-V^ M (i/) = 0. (73) 

They will be useful tools below. 

Step 2. 

If we shrink bs to b^ n S uniformly in time, transport and Gauss the- 
orems (see also (14611 and (|47|l ) allow us to obtain the pointwise balance 

-[Q]U + [S] • m + Div E X = 0, (74) 

thanks to the arbitrariness of bs- Of course, the sole difference between (|74ll 
and (|48[) is the term Div^X accounting for the interfacial structure of the 
surface S. 

Step 3. Deduction of the referential interfacial balance of standard in- 
teractions Ijfifljl . 

If f 2 acts alone, then 

X = T T v, Q = p±-v, 3 = -P T v. (75) 

Moreover, thanks to (|72(l we get 

Div s X = v-Div s T. (76) 

Then, from (|74|l we obtain Q63|l thanks to the arbitrariness of v, which is 
continuous across E. 

Step 4- Deduction of the referential interfacial balance of substructural 
interactions (|64|l . 

If G acts alone, then 

X = S T £ M (v), Q = p di,x< M W), $=-S T i M {v). (77) 

Moreover, thanks to (|73|l we get 

Div s X = £ M (i/) • (DivsS - 3) . (78) 

Then, from l|74|l we obtain (j64(l thanks to the arbitrariness of the element 
£ selected in the Lie algebra of G. 

Step 5. Deduction of the balance of configurational forces along the nor- 
mal m {i.e. (Ip3in . 

If f 1 acts alone, then 

Q = - P F T ± ■ w-p (Vi/)* Box • w, (79) 
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ff = ( Qp |x| 2 + p oX (v, *)) I-p) w (80) 

and, after some algebra, 

X= -C^W-CWm, (81) 

with Ctan and c denned respectively by (|66|l and Q67JI and w m = w • m. 

Now, equation (|74|) comes into play: we will evaluate the component 
along m of the right-hand side term of (|74|l . a vector, by taking also into 
account the arbitrariness of w. 

First we write 

~[Q]U + m ■ m = Po [F T ±]U • w+p [(Viv)* dox]U • w+ 

1 2 

+ 2/ 9 o[l*l ]w • m + p [x (i/, i>)]w • m - [P]w • m. (82) 

Moreover, by using (|7l"1) and Lemma 1, we also get 

Div s (Cf an w-(-cw m ) = w- (Diw^Ctan + (Div s c) m) . (83) 

Of course, in obtaining JSSJl, properties 3 and 4 of the definition of the 
relabeling f l of £>o including S play a crucial role. 

By inserting (|82|l and l|83|) in (|74|l . thanks to the arbitrariness of w, we 
obtain 

Po [F T ±]U + Po [(Vv)* d 0X ]U + ip [|x| 2 ]m+ 

+ Po [ X (u, v)]m = [P T ]m + DivzCtan + (Div S c) m (84) 

and we shall evaluate the component along m of (|84|l . 

First we focus our attention on terms involving x and i>. Let us introduce 
the averaged velocity v given by 

v = (4) + U (F> m. (85) 

With the help of the relation [x] = — J7[F]m introduced previously, we then 
get 

Po [F T x]t/.m = p [x].v-p [|x| 2 ]; (86) 

in other words the normal component of the vector p [F T x]C7 is equal to 
minus the jump of the relative kinetic energy \\p |x — v| ] as it is simple 
to verify. 

Still taking into account the relation [x] = — J7[F]m and the definition 
of v, we also find 

i Po [|x| 2 ] = - Po [x].v+if/ 2 [|Fm| 2 ]. (87) 

Now, by evaluating the normal component of iJSJl, using (JSSJ, (|S7Jl and 
taking into account that m-Divs < Ctnn — Ctan • L as a consequence of Lemma 
2, we obtain Ij65|l and the theorem is proven. 
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